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2.  Current  and  Anticipated  Results 

The  ultimate  objective  of  our  research  is  to  provide  algorithas,  and, 
by  using  experimental  digital  computer  programs  which  iapleaent  the  con- 
straint aethod  for  three  diaensional  analysis,  to  demonstrate  the  effec- 
tiveness of  the  method  in  three  diaensions.  Since  the  approaches  to  be 
used  in  deriving  these  algorithas  are  based  on  similar  approaches  for  two 
diaensional  analyses,  part  of  the  grant  period  has  been  spent  in  complet- 
ing the  formulation  of  algorithms  in  two  dimensions  so  as  to  have  a firm 
foundation  for  the  more  complex  three  diaensional  work.  Also  in  response 
to  a suggestion  by  Dr.  V.  B.  Venkavva  of  the  Air  Force  Flight  Dynamics 
Laboratory  at  Wright-Patterson  AF3,  the  plate  bending  element  of  the  Con- 
straint Method  is  being  developed  further.  Results  in  the  first  three 
areas  described  below  are  for  two  dimensional  problems  in  linear  stress 
analysis;  results  in  the  fourth  area  are  for  three  diaensional  linear 
stress  analysis. 

’.1  C®  Disolacement  Fields  (olane  elasticity 


.1.1  Hierarchical  Elements  and  Precomputed  Arravs 


It  has  been  shown  previously  by  others  that  elemental  arrays  may  be 
efficiently  generated  through  the  use  of  "precomputed"  arrays  - - that  is 
arrays  which  are  computed  once,  stored  on  permanent  file,  and  then  reused 
in  all  subsequent  applications  of  the  program.  Hie  new  work  done  by  the 
principal  investigator  and  his  collaborators  has  two  objectives:  the  first 
is  to  show  how  the  hierarchal  C*  elements  (described  in  1.2)  for  a quadratic 
functional  aay  be  formulated  using  precomputed  arrays  thus  yielding  a fi- 
nite element  technique  which  is  especially  suited  to  problems  with  local 
rapid  variation  of  the  function  to  be  approximated.  In  particular,  formulas 
for  two- diaensional  (hierarchal)  elenenc  arrays  for  arbitrary  polynomial 
order  are  derived,  based  on  precomputed  arrays.  The  second  objective  is  to 
apply  the  combined  approach  of  hierarchal  elements  and  precomputed  arrays  to 
decide  if  a computed  result  has  "converged".  A common  practice  in  finite 
element  analysis  is  to  solve  a problem  several  times  using  successively 
required  meshes  i.e.  to  apply  the  procedure  for  h-convergence . If  suc- 
cessive analyses  agree  then  it  is  usually  assumed  that  the  finite  element 
approximation  is  accurate.  This  procedure  can  be  computationally  expensive 
when  several  highly  refined  meshes  are  used.  An  alternative  procedure  is 
to  use  p-convergence,  which,  as  pointed  out  in  1.1,  has  a faster  rate  of 


convergence  to  the  true  diaplecemencs . The  computational  effectiveness  of 
the  p-convergence  procedure  is  demonstrated  numerically  using  hierarchal 
elements  and  precomputed  arrays. 

Detailed  formulas  are  given  for  calculation  of  stiffness  matrices  and 
for  calculation  of  polynomial  coefficients  from  nodal  variables,  Hierar- 
chal  nodal  variables  are  presented  together  with  some  of  the  favorable  con- 
sequences of  using  hierarchal  nodal  variables.  Computation  times  for 
stiffness  matrices  are  given  in  terms  of  equivalent  time  units  (e.t.u.)  for 
different  methods.  It  has  been  demonstrated  in  [Li]  that  the  constraint 
method  allows  many  more  degrees  of  freedom  than  other  mechods  do.  for  the 
same  computer  cost. 

2.1.2  Singularity  Functions  used  In  Linear  Elastic  Fracture  Mechanics 

Poor  computational  efficiency  has  generally  been  observed  when  employ- 
ing conventional  finite  elements  near  a crack  top.  Therefore,  in  [12],  the 
behavior  of  a rational  type  singularity  function  which  represents  an  order 
of  r 3 tvpe  stress  singularity  has  been  investigated,  when  used  in  con- 
junction with  hierarchical  C*  elements,  in  an  effort  to  improve  erficiencv. 
An  internal  mode  of  the  form 


Li  L:  L3 

*ai*  L:’  L3}  ■ (L,  + T.J J/ 2 
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is  an  approximation  in  triangular  coordinates  to  an  r * type  stress  singu- 
larity at  the  vertix  ■ 1.  Using  techniques  similar  to  those  developed  by 
us  in  [5],  it  is  possible  to  Integrate  derivatives  of  the  singularity  func- 
tion J over  triangles,  explicitly,  and  then  to  employ  precomputed  arrays  to 
compute  elemental  stiffness  macrices.  The  p-convergenc  procedure  can  then 
be  applied  to  determine  the  strain  energy  U,  the  strain  energy  release  race 
C and  crack  opening  displacements. 

To  determine  the  effect  of  adding  the  singularity  function  » to  a 
polynomial  basis  we  heve  taken  the  edge  cracked  panel  shown  in  figure  6. 
Because  of  symmetry  only  one  quarter  of  the  panel  is  modeled  with  Che  three 
different  triongulatlons  shown  in  figure  6.  The  finite  element  solution 
is  obtained  for  each  of  the  triangulations  firsc  with  only  polynomials  in 
the  displacement  field,  and  then  with  the  addition  of  the  singularity  func- 
tion in  elements  meeting  at  the  crack  cop.  The  results  ars  plotted  In 


figure  6 and  figure  7 and  are  given  in  decail  in  Table  2 and  Table  3.  Ic 
la  evident  char  there  is  an  laprovenent  in  the  estimates  of  the  strain 
energy  U and  the  strain  energy  release  rate  G but  that  this  laproveaent  is 
not  appreciable  coopered  to  the  results  already  obtained  by  increasing  the 
polynomial  order  p.  Thus,  the  mode  of  p-convergence  appears  to  be  more  im- 
portant in  achieving  accuracy  chan  the  addition  of  a singularity  function 
to  the  basis. 

Ve  have  also  considered  the  case  of  a centrally  cracked  panel  shown 
In  figure  3 aodeled  with  five  eleaents  in  one  quarter  of  Che  panel.  The 
crack  opening  displacement  (CCD)  and  stresses  are  computed  with  the  singu- 
larity function  included  in  the  approximating  displacement  field  and  the 
results  are  shown  in  table  4.  In  Cable  4,  5 is  Che  CCD  obtained  by  polyno- 
mial approximation  above,  5*  is  che  COD  obtained  by  including  a singu- 
larity function  and  c la  a measure  of  component  defined  as 

* 

c - x 100  . 

It  is  clear  froo  table  4 that  inclusion  of  the  singularity  function  in  the 
p-esnvergent  procedure  leads  to  substantial  improvement  in  che  displacement 
near  the  crack. 

More  details  and  ocher  cases  which  have  been  studied  are  presented  in 
[17],  In  all  cases  che  case  of  the  constraint  method  in  the  mode  of  ? 
convergence  leads  to  major  computational  advantages. 

2.2  Co'ioled  and  C Displacement  Fields 

la  this  work  Che  results  of  Xratochvii  et  al  in  [13]  are  generaliced 
to  problems  with  three  independent  displacement  fields,  .in  essential  aspect 
of  this  approach  is  to  transform  a triangular  element  T in  the  x-y  plane 
into  a standard  triangle  t with  vertices  at  the  origin  and  at  a unit  dis- 
tance along  che  horizontal  and  vertical  axes.  Such  a transformation  is 
shown  on  che  next  page  where  I and  i represent  coordinates  in  che  plane  of 
the  standard  triangle  and  x and  y are  local  coordinates  for  che  elemenc 
with  the  first  vertex  of  the  element  coinciding  with  the  origin.  The  other 
two  vertices  and  also  che  three  edges  are  numbered  in  counter-clockwise 
order  as  shown. 


The  justification  for  using  the  standard  triangle  is  that  integrations  and 
matrix  inversions  are  performed  with  respect  to  the  standard  triangle.  Thus 
they  need  be  dene  only  once  and  the  results  are  stored  and  then  used  in 
all  future  applications  of  the  program.  Computation  of  the  element  stiffness 
aatrix  is  thus  reduced  to  computing  a linear  combination  of  a small  number 
of  precomputed  matrices  followed  by  pre-  and  post-multiplication  by  block 
diagonal  matrices.  The  number  of  precomputed  matrices  which  must  be  stored 
are  considerably  reduced  by  choosing  hierarchal  nodal  variables.  Further 
details  are  given  in  [12]. 

2.3  C~  Displacement  Fields  'dace  bending  problems) 

It  was  suggested  by  Dr.  V.  3.  Venkayya  (of  the  Analysis  and  Opcimi- 
ration  Group.  Structures  Division,  Air  Force  Flight  Dynamics  Laboratory, 
Vright-?attar3on  Air  Force  3ase)  in  a letter  dated  12  December  1975  that 
further  development  of  place  bending  elements  would  be  useful  in  order  to 
realize  the  full  potential  of  the  constraint  method.  Accordingly,  a 
sophisticated  place  bending  elemenc,  incorporating  a complete  order 
polynomial  with  ? i 5 and  corrective  rational  functions,  has  been  formu- 
lated and  is  now  being  programed  and  cescad.  Ve  now  describe  some  of 
this  work. 

It  iJ  well  known  (see  [19],  for  example)  that  exactly  conforming 
(even  at  vertices)  C1  displacement  fields  cannot  be  formed  merely  by  freely 
assembling  finite  elements.  There  are  certain  additional  constraint  equa- 
tions which  must  be  satisfied  at  vertices.  The  simplest  form  of  these 
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constraint  equations  has  been  given  by  Peano  in  [4],  where  a specially 
devised  assembly  procedure  is  presenced  which  automatically  enforces  the 
constraint  equations.  An  alternative  method  for  enforcing  exact  conformity 
when  using  displacement  fields  of  arbitrary  polynomial  order  p is  by  sup- 
plementing p^-  order  polynomials  with  newly  constructed  corrective  rational 
tunctions.  This  destroys  the  analytic  character  of  the  approximation  at 
the  vertices  but  permits  free  assembly  of  elements  without  enforcing  con- 
straints [4],  An  algorithm  has  been  developed  for  a (exactly)  conform- 
ing triangular  element  which  contains  complete  polynomials  of  order  p > 5 
and  corrective  rational  functions  A typical  element  of  order  p * 5 to- 
gether with  nodal  variables  is  shown  below. 


3 


The  quintic  element  has  24  independent  nodal  variables.  The  shape 

functions  for  these  nodal  variables  are  given  in  Table  5.  The  shape  func- 
tions for  second  order  tangential-normal  derivatives  are  rational  functions. 
It  is  important  to  observe  that  although  rational  functions  are  used  in 
the  basis,  all  terms  which  appear  in  the  elemental  stiffness  matrix  can  be 
integrated  explicitly  without  recourse  to  numerical  quadrature.  This  was 
proved  in  [5].  An  algorithm  based  on  a hierarchical  family  of  elements 
using  corrective  rational  functions  has  been  programmed  and  is  now  being 
tested  on  numerical  examples.  We  now  give  one  such  example. 


2.3.1  An  Example:  A Simolv  Supported  Square  Plate  under  Uniform  Loading 

Results  obtained  by  using  the  hierarchical  quintic  C1  element  are 
compared  with  those  obtained  by  Covrper  in  [20],  by  Caramanlian  et  al  in  [21] 
and  by  Tsai  in  [22].  Two  different  types  of  triangulation  are  shown  in 
Figure  9:  the  Q-arrangement  and  the  P- arrangement.  In  Table  6 the 
results  for  central  deflection,  central  bending  noments  and  for  the  strain 
energy  are  compared  for  different  elemental  arrangements  and  for  different 
elements.  It  is  seen  that  for  the  Q-arrangement  the  strain  energy  ob- 
tained by  using  the  Constraint  Method  is  an  order  of  magnitude  more  accurate 
than  that  obtained  by  Cowper,  and  for  the  p-arr3ngement  it  is  two  orders 
of  magnitude  more  accurate.  This  is  consistent  with  the  order  of  magnitude 
improvements  that  have  resulted  in  many  problems  by  using  the  constraint 
method. 

2.4.  Hierarchical  Families  of  Complete  Conforming  Solid  21ements  of 
Various  Shapes  and  Arbitrary  Order 

In  [4]  a table  of  canonical  basis  functions  for  a triangular  (two 
dimensional)  element  was  presented,  using  natural  coordinates.  This  table 
has  been  ger.eraliced  to  include  a canonical  basis  for  a tetrahedral  (three 
dimensional)  hierarchic  family.  Using  this  table  shape  functions  are  gene- 
rated for  tetrahedral  C°  elements  and  their  corresponding  nodal  variables. 
Table  7 gives  the  nodal  variables  and  shape  functions  for  the  first  four 
hierarchic  Tetrahedral  C®  elements. 

A hierarchic  family  of  rectangular  C®  elements  has  been  developed  and 
this  family  is  used  to  generate  a hierarchic  family  of  C®  brick  elements 
of  arbitrary  polynomial  order.  3y  using  a combination  of  the  hierarchic 
triangular  and  rectangular  families,  we  have  also  constructed  a hierarchic 
family  of  triangular  prismatic  elements.  These  prismatic  elements,  which 
have  also  have  arbitrary  polynomial  order,  can  be  made  to  join  continuously 
to  tetrahedral  elements,  so  that  pointed  prismatic  geometries  can  now  be 
easily  approximated  (see  the  figure  on  the  next  page) . 
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TABLE  2 

Comparison  of  Strain  Energy  Approximations 
for  Various  Polynomial  Orders  and  Mesh  Divisions 
With  and  Without  the  Inclusion  of  Singularity  Mode 


p 

Mesh 

•)  n 

U E/a”l“ 

P 

Strain 

Error 

(S) 

Energy 

U*E/o“l 

P 

Error 

(*) 

A 

2.9532796 

3.356 

2.9815737 

7 . 634 

2 

3 

2.9773692 

7.749 

2.9909923 

7.342 

C 

2 . 9876714 

7.445 

2.9957397 

7.194 

A 

3.0809992 

4.554 

3.0976153 

4.039 

3 

3 

3.0961462 

4.035 

3.1050224 

3.310 

C 

3.1124021 

3.531 

3.1172042 

3.432 

A 

3.1391130 

2.753 

3.1503637 

2.405 

4 

3 

3.1433200 

2.453 

3.1547416 

2.269 

C 

3.1603298 

2.096 

3.1635274 

1.997 

A 

3.1683329 

1.347 

3.1761374 

1.607 

5 

3 

3.1751263 

1.638 

3.1791517 

1.513 

C 

3.1326309 

1.405 

3.1343030 

1.333 

A 

3.1350973 

1.329 

3.1907447 

1.154 

6 

3 

3.1399564 

1.179 

3.1923865 

1.  OSS 

C 

3.1955412 

1.001 

3.1971365 

0.956 

♦Singularity  function  is  included  in  the  elements 
meeting  at  the  crack  tip. 


-17- 


TAELE  3 

Comparison  of  Strain  Energy  Release  Race  Approximations 
for  Various  Polynomial  Orders  and  Mesh  Divisions 
With  and  Without  Che  Inclusion  of  Singularity  Mode 


Strain 

Energy 

Release  Rate 

p 

Mesh 

G E/c“i 

P 

Error 

(2) 

G*E/o“l 

P 

Error 

C) 

A 

2.2735340 

13.510 

2.3195304 

16.362 

2 

B 

2 . 3543444 

15.597 

2.3663902 

15.165 

C 

2.3926472 

14.242 

2.4012298 

13.934 

A 

2.5083682 

10.094 

2.537S67S 

9.043 

3 

B 

2.5660932 

3.025 

2.5726530 

7.790 

C 

2.5995090 

6.323 

2.6035322 

6.632 

A 

2.6206222 

6.071 

2.6405090 

5.353 

4 

3 

2.6592740 

4.636 

2.6636310 

4.529 

C 

2.6330336 

3.334 

2.6859452 

3.729 

A 

2.6752748 

4.112 

2.6883926 

3.642 

5 

3 

2.7015356 

3.169 

2.7041712 

3.076 

C 

2.7134304 

2.565 

2.7203946 

2.495 

A 

2.7074112 

2.960 

2.7170000 

2.616 

6 

3 

2.72625*0 

2.235 

2.7232633 

2.213 

C 

2.7391374 

1.823 

2.7405716 

1.772 

"Singularity  function  is  included  in  the  elements 
meeting  at  Che  crack  tip 
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TA3L2  4 


Comparison  of  tha  Crack  Opening  Displacements  (COD) 
for  Various  Polynomial  Orders 
With  and  Without  the  Inclusion  of  Singularity  Modes 


«/«.  :S*  llatasca  ilMf  m crack  wriic)  !roo  ::<  caocar  a<  craca 

COB 


9.0 

O.i 

O.I 

0.3 

9.1 

0.3 

0.4 

0.7 

O.i 

0.9 

s 

9..1M 

3.. 440 

3.-434 

3.- 419 

9.13*1 

o..;:3 

0.3377 

:.:ss9 

0.1030 

0.1100 

3.30*2 

3.3300 

3.51i: 

3.3074 

0.-791 

3.-342 

0.37)1 

0.3039 

0.1137 

0.1147 

1 

1.5 

1.3 

5.1 

3.3 

3.3 

3.7 

3.3 

3.9 

« a 

4.1 

i 

3.43:* 

9.9304 

0.3071 

3.3404 

3.343: 

3.<:44 

0.-714 

0.3933 

o.i?*: 

0.1454 

0.4—3 

9.4304 

0.4171 

0.40*0 

0.341- 

0.3*90 

0..424 

0.-143 

0.3111 

0.1750 

< 

l.t 

1.6 

1.4 

• • • 

1.9 

3.4 

a(a 

3.3 

4.1 

4.3 

i 

3.479) 

3.449* 

0.43:: 

0.430* 

3.3034 

3.5440 

3.3:1: 

o.iiia 

0.3317 

0.1041 

3.9445 

3. *7*1 

3.3349 

0.4343 

3.4101 

0.3779 

3.3337 

3. -443 

0.3449 

3.1193 

C 

i.O 

1.3 

1.3 

l.l 

1. 6 

a ^ 

2 . J 

*.i 

i 

3.4913 

3.94-3 

3.i7:r 

3.45:7 

3.41:7 

3.344! 

0.3397 

3..  77: 

3.3431 

3.5134 

)• 

3.4439 

0.4449 

3.47)* 

3.9541 

3.9131 

0.590* 

0.5-51 

3.-  443 

0.5943 

O.DU 

< 

3.7 

3.7 

O.i 

0.4 

0.7 

0.7 

1.3 

1.9 

3.3 

3.4 

i 

0.447) 

3.49t3 

3 . 44C0 

0.3433 

0.9339 

3.5943 

3.5*40 

0.-444 

0 • *016 

3.1421 

i • 

3. 7313 

3.4454 

3.  :43* 

0. 900- 

0. 

3.9014 

3.53:0 

3.-911 

3.-133 

0.17— 

€ 
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0.3 

3.9 
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3.4 

3.3 

3.4 
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-.7 

i 
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3.4443 

3. 4444 

3.4447 

3.4-03 

3.40*4 

3.3354 

3. .91: 

3. -391 

o.i:): 

i* 

3.7039 

3. '31* 

3.9493 

0.4441 

3.4*3* 

3. 408* 

3.5390 

0. -4*5 

3.-159 

3.133* 
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0.4 

3.4 

0.3 
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TABLE  D 

SHAPE  FUNCTIONS  FOS  FIFTH  onDE.'t  C1  ELEMENT 

CONTAINING  CORRECTIVE  RATIONAL  FUNCTIONS 


3 > 2 

l2-l2 

m 1 3 

‘ 


NODAL  VARIABLE 

-(3’  ) 


3Sw 


3*33b3 
2 


■a> 


»\ 


,2 

j v 


’V"2 


a) 


SHAPE  FUNCTION 


*21 


.2  .2 

1 “2  “1 

1 .2  .2 


S4  " J “1  “2 

N,  ■ “ 1?  L3U-L,)  - 

O * 1 J J «*  * 

♦ r<i-u  }?!„ 

* • * 

(N,fl  - 1?  L2  tj) 


* A J 


r— <1) 
J33 

m 

SJ2 

w(l) 


3*33n3 


(1) 


3a. 


‘S-ui 


i,n3 


»2  - iJu-tj).a5-«»4 

"j  - ij  ljU-1,1 

N,  - L3  *3(N,*N,)-6(M  «N„-N 
1 1 2 3 4 3 

L‘  1,3  r 

1 . II -3 

1 L,  * L, 


1 . £ h •’ 

J L,  - 


MJ 


CORRECTIVE 

RATIONAL 

FUNCTIONS 


CENTRAL  UKNIHNC 

MOMENTS  CENTRA  I.  DEFLECTION 
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TABLE  6 

SIMPLY  SUPPORTED  SQUARE  PLATE 
UNIFORM  LOAD 

(5ctl  order  Approximation) 


N 

P/Q 

Constraint 

Method 

Covper 

Caramanlian 

Tsai 

* 

Q 

4.060210 

4.0609374 

4.060236 

4.0597 

? 

4.059696 

4.0634349 

4.069917 

4 

Q 

4.062342 

4.0623473 

4.062323 

4.0623 

P 

4.062433 

4.0627265 

4.06250 

6 

Q 

4.0623493 

4.0623517 

P 

4.0623503 

4.0623393 

3 

Q 

4.0623522 

4.0623524 

EXACT  SOLUTION 

: 4.0623527 

N 

P/Q 

Constraint 

Method 

Covper 

Caramanlian 

Tsai 

2 

Q 

4.63094 

4.70261 

4.6579 

4.6127 

p 

4.34074 

4.31530 

4 . 3563 

4 

Q 

4.73697 

4.73259 

4.7316 

4.3333 

p 

4 . 790  74 

4.79033 

4.7924 

6 

Q 

4.737436 

4.737621 

P 

4.739023 

4.738972 

3 

Q 

4.738329 

4.738419 

EXACT  SOLUTION:  4.7'3633 


. 
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TA3LE  6 (continued) 


STRAIN  ENERGY 
104  TO/qQL6 


N 

P/Q 

Cowper 

" Error 

Constraint 

Method 

Z Error 

2 

Q 

3.5099612 

3.04x10“ 2 

3.5117139 

9 . 79*10*" 3 

P 

3.4776356 

4 .10*10~2 

3.5124319 

3.31xl0“4 

4 

Q 

3.5124403 

1. 32*10~3 

3.5125393 

1.56xl0“4 

P 

3.5113962 

1.36xl0-2 

3.5125227 

3.51xl0“4 

6 

Q 

3.5125336 

1.64xl0~4 

3.5125514 

1.41xl0“5 

P 

3.5124190 

1.57xl0“3 

8.5125501 

2.94xl0~5 

3 

Q 

3.5125496 

3.52xl0“3 

3.5125524 

2.35xl0“6 

EXACTION  SOLUTION:  3.5125526 


hOOAl  VAIUAOIES  AMD  SHAPE  EUMCTIOMS  FOR  THE 
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(D.Sc  May  1978,  Thesis  Title:  "p-convergent  finite 
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INTERACTIONS 


On  October  19,  1977,  the  Principal  Investigator,  I. 
Norman  Katz,  presented  an  invited  talk  at  the  Expository 
Seminar  Series  of  the  Applied  Mathematics  Division  of  the 
U.S.  National  Bureau  of  Standards  in  Gaithersburg,  MD.  An 
abstract  of  the  talk,  based  on  current  research,  is  enclosed. 


Mac i or. a 1 3ureau  of  Standards 
Applied  Mathematics  Division 

Exnosicorv  Seminar  Series 


0CT03.r3.  meeting 


Data  & Time:  Wednesday , October  19,  1977 11:00  A.M. 

(Coffee  Social  10:45  A.M.) 

Place : Technology  Building,  Room  327 

National  3ureau  of  Standards 
Caichersburg,  Maryland  20760 

Speaker:  Professor  I.  Norman  Fate 

Department  of  Systems  Science  and  Mathematics 
Washington  University 
Sc.  Louis,  Missouri 


Title  of  Talk: 


"The  Constraint  Method  tor  Finite  Element  Stress  Analysis" 


Abs trac t : In  conventional  approaches  to  finite  element  stress  analysis 

accuracy  is  obtained  by  fixing  the  degree  p of  the  approximat- 
ing polynomial  and  by  allowing  the  maximum  diamecer  h of  ele- 
ments in  the  triangulation  to  approach  rsro.  An  alternate 
. approach  is  to  fix  the  triangulacion  and  to  increase  the  degrees 

of  approximating  polynomials  in  those  elements  where  more  accuracy 
is  required.  In  order  to  implement  the  second  approach  efficiently 
it  is  necessary  to  have  a family  of  finite  elements  of  arbitrary 
bolyr.om.ial  degree  p with  the  property  that  as  much  information 
as  possible  can  be  retained  from  the  pm*I  degree  approximation  when 
computing  the  (p-i-l)sc  degree  approximation.  Such  a HlZILARC'dlC 
family  has  been  formulated  with  p > 2 for  problems  in  plane 
scress  analysis  and  with  p > 5 for  problems  in  plate  bending. 

The  family  is  described  and  numerical  examples  are  presented  which 
illustrate  the  efficiency  of  the  new  method. 

3io  wrmhical 

ike ton  : Professor  'Lit:  received  his  Ph.D.  in  Mathematics  from  M. 1.7.  in 

1959.  "He  has  worked  at  AVCO/Research  and  Advanced  Development 
in  Wilmington,  Massachusetts  where  he  became  Manager  of  the  Mathe- 
matics Department.  Since  1967  he  has  been  at  Washington  University 
in  St.  Louis.  Ris  research  has  been  in  numerical  analysis,  ordinal- 
and  partial  differencial  equations,  finite  elemenc  mechods,  optimal 
facility  locacion  and  biocachemacics . 


All  inti  res  ted 


a re 


invited  to  actend. 


